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I.  Introduction 


The  objectives  of  our  project  funded  by  BMDO/IST  as  managed  by  ONR  under 
the  direction  of  Dr.  Richard  S.  Miller  were  the  following: 

1.  To  develop  a  source  code  to  mathematically  model  the  thermo-mechanical 
history  experienced  by  energetic  materials  during  continuous  processing  in  a  co¬ 
rotating  twin  screw  extruder. 

2.  To  experimentally  validate  the  mathematical  model  using  well-instrumented 
industrial-scale  twin  screw  extruders. 

The  project  satisfied  these  objectives  as  well  as  discovering  the  following  which 
have  significant  impact  on  the  continuous  processability  of  energetics. 

•  A  particle-free  region  develops  adjacent  to  all  walls  during  the  flow  of  energetic 
materials  with  their  characteristic  high  degree  of  fill  with  solids. 

•  This  particle  free  and  binder-risk  region  at  the  wall  gives  rise  to  a  region  with  a 
high  -deformation  rate  and  thus  serves  as  an  apparent  slip  layer. 

•  This  thickness  of  this  layer  can  be  determined  along  with  the  wall  slip  velocity 
versus  the  imposed  shear  stress  condition  by  using  viscometric  flows. 

•  The  wall  slip  behavior  needs  to  be  incorporated  into  mathematical  models  of  the 
process. 

•  The  formation  of  this  wall  slip  layer  is  affected  by  the  presence  of  entrained  air 
which  further  lubricates  the  surface. 

•  The  binder  polymer  and  the  solids  can  separate  along  the  flow  direction  where 
the  basis  of  the  demixing  is  the  filtration  of  the  binder  in  the  flow  direction.  This 
filtration  phenomenon  is  a  well-defined  safety  hazard  in  processing  of  energetics. 

•  The  amount  of  air  entrained  during  processing  is  affected  by  screw  geometries 
and  operating  conditions  used  in  extrusion. 


•  Extruders  can  develop  hot  spots.  With  increasing  degree  of  fill  of  the  energetic 
suspension,  the  relaxation  time  of  the  suspension  increases  (it  becomes  more 
solid-like)  resulting  in  a  yield  stress  (viscoplasticity),  which  exasparates  the 
formation  of  the  hot  spot.  The  development  of  the  hot  spot  is  a  safety  hazard 
which  should  be  eliminated. 

•  Twin  screw  extruders  can  convey  energetic  materials  without  mixing  them 
giving  rise  to  "pipeline"  flows.  The  formation  of  such  flow  channels  for  unmixed 
materials  constitutes  a  safety  hazard  and  deteriorates  the  mixing  efficiency  of  an 
extruder. 

Our  discovery  of  these  phenomena  specific  to  the  highly  filled  nature  of  energetic 
materials  lead  us  to  develop  tools  (both  experimental  and  theoretical)  to  elucidate 
them  further.  In  the  following,  the  understanding  which  we  have  gained  from  our 
studies  of  the  flow  and  deformation,  continuous  processability,  modeling  and 
microstructure  development  aspects  of  processing  of  energetic  suspensions  are 
summarized.  The  publications  which  resulted  from  each  area  of  study  are  also 
included  and  they  can  be  consulted  for  further  information.  Finally,  the 
processing  configurations  "lessons  learned"  from  our  studies  are  also  included. 

I.I.  Issues 

Manufacture  of  propellants,  explosives  and  pyrotechniques  (PEP's)  necessitate 
the  processing  of  a  energetic  suspension  or  a  dispersion  at  a  solid  concentration 
which  approaches  the  maximum  packing  fraction,  <t)m,  of  the  solid  phase.  The 
concept  of  the  maximum  packing  fraction  is  schematically  illustrated  in  Figures  1- 
3.  The  use  of  a  multimodal  size  distribution  (Figure  2)  is  shown  to  increase  the 
amount  of  solids  which  can  be  incorporated  into  a  given  volume  i.e.,  a  higher 
maximum  packing  fraction  in  comparison  to  a  unimodal  size  (Figure  1),  whereas 
the  use  of  non-spherical  and  fibrous  particles  (Figure  3)  decreases  the  maximum 
packing  fraction  [1]. 

The  use  of  particles  with  different  sizes  is  necessary  in  formulation  of  energetic 
materials  to  increase  the  maximum  packing  fraction.  This  is  illustrated  further  in 
Figure  4.  Here  the  maximum  packing  fraction  is  shown  to  be  affected  by  the  size 
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Figure  4.  Maximum  packing  fraction  versus  size 
concentrations  of  spherical  particles  at  different  size 
ratios,  a,  for  bimodal  distribution:  experiments  and  theory. 


ratio,  a,  i.e.,  ratio  of  diameters  of  the  large  over  small  spheres,  and  the 
percentage  of  the  smaller  spheres  in  a  bimodal  particle  size  distribution  [2,  3]. 
Thus,  the  shape  and  size  distributions  of  solid  particles  are  critical  in  determining 
the  maximum  possible  concentration  of  solids  incorporated  into  a  suspension 
consisting  of  a  liquid  matrix  phase  (the  binder)  and  solid  particles. 

The  common  denominator  in  propellant  formulations  is  the  use  of  the  minimum 
amount  of  the  liquid  phase  "binder"  necessary  to  facilitate  the  flowability  and 
shapeability  of  the  energetic  suspension  and  the  maximum  amount  of  solids  to 
achieve  the  best-possible  combustion  properties.  The  minimal  use  of  the  binder 
by  working  close  to  the  maximum  packing  fraction  of  the  energetic  solid  phase 
reduces  the  porosity  in  the  final  grain.  Overall,  the  high  degree  of  fill  of  such 
highly  concentrated  energetic  suspensions  introduces  significant  difficulties  in 
their  characterization,  processing,  simulation,  and  quality  control.  However, 
propellants,  regardless  of  their  application,  exhibit  a  number  of  common  flow, 
deformation  and  processing  traits  which  can  be  understood  and  used  for  their 
optimum  design  and  operation  of  continuous  manufacturing  operations.  In  this 
final  report,  various  hitherto  identified  traits  will  be  outlined  and  ramifications  on 
processing  operations  will  be  provided.  The  report  will  also  review  the  simulation 
tools  which  we  have  developed. 

As  shown  schematically  in  Figure  5,  our  research  group  has  developed 
capabilities  in  simulation  (using  three  dimensional  Finite  Element  Method  based 
source  codes  developed  in  house),  rheological  characterization,  microstructural 
analysis,  wettability  and  work  of  adhesion  between  solid  and  liquid  phases,  and 
characterization  of  ultimate  properties  of  propellant  formulations  and  their 
simulants.  In  the  following,  some  of  these  experimental  and  numerical  analysis 
tools  will  be  described,  along  with  some  of  the  technologies  developed. 


III.  Flow  and  Deformation  Behavior  of  Propellants  and  their  Simulants 


The  characterization  of  the  flow  and  deformation  behavior  of  propellants  which 
are  highly  filled  suspensions  is  first  and  foremost  complicated  by  the  occurrence 
of  wall  slip  [4-16].  The  wall  slip  behavior  of  a  concentrated  suspension  consisting 


Figure  5.  The  integration  of  numerous  capabilities  developed. 


of  sixty  three  percent  by  volume  of  glass  spheres  in  an  acrylonitrile  terminated 
polybutadiene,  PBAN,  matrix  is  demonstrated  in  Figure  6.  The  suspension 
sample  undergoes  steady  torsional  flow  in  between  two  parallel  disks;  the  disk  at 
the  top  is  rotating  and  the  disk  at  the  bottom  is  stationary.  Before  the  onset  of 
deformation,  a  straight  line  marker  is  placed  at  the  free  surface  of  the  suspension 
and  the  edges  of  the  two  disks.  Upon  the  rotation  of  the  top  disk,  discontinuities 
appear  at  both  the  top  and  bottom  suspension/wall  interfaces,  indicating  wall  slip. 

This  technique  which  we  have  developed  allows  the  unambiguous  determination 
of  the  slip  velocity  and  the  true  deformation  rate  of  the  suspension  directly  [13]. 
Computerized  image  analysis  of  the  recorded  images  generates  the  velocity  of 
the  suspension  found  adjacent  to  the  wall,  V,  as  a  function  of  time.  The  slip 
velocity  at  the  wall,  0S,  is  0S  =  Vf  -Vw,  where  Vw  is  the  wall  velocity.  This  is 
one  of  the  techniques  we  have  used  to  characterize  the  wall  slip  behavior  of 
propellants  and  their  simulants. 

In  Figure  7  the  time  dependent  development  of  the  wall  slip  behavior  of  a 
suspension  in  steady  torsional  flow  is  shown  [16].  With  time,  the  deformation  rate 
of  the  suspension  decreases,  while  the  shear  stress  and  the  wall  slip  velocity 
increase.  When  the  steady  state  behavior  is  achieved,  the  true  deformation  rate 
is  only  about  1/5th  of  the  imposed  apparent  shear  rate.  The  thickness  of  the  slip 
layer  can  be  estimated  as  a  function  of  the  viscosity  of  pure  polymer  binder  and 
is  generally  in  the  order  of  microns  [9,  13,  16]. 

If  the  surface  to  volume  ratio  in  rheometry  is  systematically  changed,  the  slip 
velocity  can  be  determined  also  from  the  slope  of  the  apparent  shear  rate  versus 
reciprocal  diameter  or  gap  data  collected  at  constant  wall  shear  stress.  For 
example,  in  steady  parallel  disk  torsional  flows  the  apparent  shear  rate  at  the 
edge,  7aR>  becomes  proportional  to  reciprocal  of  the  gap  between  the  two  disks, 
l/H  [7],  Plots  of  apparent  shear  rate,  taRj  versus  the  reciprocal  gap  height,  1/H, 
at  constant  values  of  shear  stress,  tr,  generate  straight  lines  with  slopes  equal  to 
2  Us(tr)  as  shown  in  Figure  8  [8,  9,  13,  14],  The  slope  hence  the  slip  velocity 
values  increase  with  increasing  wall  shear  stress.  On  the  other  hand,  in  capillary 
flow  changing  the  surface  to  volume  ratio  requires  that  the  capillary  diameter  be 
changed  while  keeping  the  length/diameter  ratio  of  the  capillary  a  constant.  Upon 
wall  slip,  the  apparent  shear  rate  varies  linearly  with  reciprocal  of  capillary 
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Figure  6.  Demonstration  of  wall  slip  of  a  concentrated 
suspension  in  steady  torsional  between  two  parallel  disks. 


Figure  7.  Time-dependent  development  of  wall  slip  in 
torsional  flow  of  a  PBAN  suspension  with  63%  by  volume 
spherical  solids. 
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Figure  8.  Apparent  shear  rate  versus  reciprocal  of  gap  height  for  a 
suspension  with  76.5%  by  volume  solids  in  steady  torsional  flow. 


diameter,  1/D,  at  constant  wall  shear  stress  with  a  slope  equal  to  8US  [4,  9], 
where  Us  is  the  wall  slip  velocity.  The  typical  steady-state  slip  velocity  at  the  wall 
versus  wall  shear  rate  behavior  of  two  propellant  simulants  are  shown  in  Figure 
9.  The  wall  slip  velocity,  Us,  varies  linearly  with  wall  shear  stress,  tr  [9,  13]. 

IV.  Solid-like  Behavior  with  Wall  Slip 

The  relative  effect  of  the  slip  flow  on  the  overall  flow  rate  in  a  die  with  an 
equivalent  diamter  of  D,  can  be  ascertained  by  the  characterization  of  the  Qs/Q 
ratio  where  Qs  is  the  volumetric  flow  rate  due  to  slip  i.e.,  Qs  =  (7t/4)D2Us  and  Q  is 
the  volumetric  flow  rate  given  by  Q  =  (7t/4)D2V.  When  the  mean  velocity  of  the 
suspension  in  the  die,  V,  becomes  equal  to  the  slip  velocity,  Us,  the  suspension 
flows  without  deformation,  that  is  as  a  plug,  lubricated  at  the  wall. 

Qs/Q  versus  wall  shear  stress  behavior  of  two  propellant  simulants  in  capillary 
flow  are  shown  in  Figure  10.  The  data  indicates  that  the  60  percent  by  vol.  solids 
propellant  simulant  flows  as  a  plug  above  a  wall  shear  stress  value  of  40  kPa, 
whereas  the  suspension  with  76.5  percent  by  vol.  solids  exhibits  plug  flow  in  the 
capillary,  at  wall  shear  stress  values  which  are  smaller  than  about  10  kPa  [9,  13]. 

These  observations  indicate  that  propellants  become  solid-like  in  various  flow 
regimes.  Solid-like  behavior  at  wall  shear  stress  values  smaller  than  a  yield 
stress  is  defined  as  viscoplasticity  and  is  related  to  the  formation  of  a  temporary 
network  which  is  destroyed  when  the  shearing  stress  exceeds  the  yield  stress. 
On  the  other  hand,  we  have  found  that  the  development  of  solid-like  behavior 
above  a  critical  wall  shear  stress  is  associated  with  dilatancy  [17]  possibly  arising 
from  the  radial  migration  of  solid  particles  away  from  the  wall  of  the  die  [18-20] 
and/or  interlocking  to  establish  solid-like  behavior. 

Figure  11  illustrates  the  viscoplasticity  of  a  propellant  simulant  with  76.5%  by 
volume  solids  [13],  This  is  the  simulant  which  was  used  to  develop  the  process 
in  the  Advanced  Shuttle  Rocket  Motor  Project  run  by  Lockheed/Aerojet 
consortium.  Upon  rotation  of  the  disk  at  the  top,  the  propellant  simulant  initially 
deforms.  However,  further  angular  displacement  of  the  rotating  disk  induces  only 
wall  slip  adjacent  to  the  walls  and  the  propellant  simulant  no  longer  deforms 
under  the  prevailing  conditions  of  constant  shear  stress.  The  inseparable  nature 
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Figure  9.  Wall  slip  velocity  versus  shear  stress  of  two  suspensions 
with  60%  and  76.5%  by  volume  solids. 


Figure  10.  Ratio  of  volumetric  flow  rate  due  to  slip,  Qs  over 
volumetric  flow  rate  for  capillary  flow,  Q,  of  two  suspensions 
with  60%  and  76.5%  by  volume  solids. 
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Figure  1 1 .  Demonstration  of  the  viscoplasticity,  i.e.,  solid-like  behavior  at  a  shear 
stress  smaller  than  yield  stress. 


of  wall  slip  and  solid-like  behavior  of  propellants  and  their  simulants  is  thus 
emphasized. 

V.  Flow  Instabilities  of  Propellants  and  their  Simulants 

We  have  discovered  that  the  pressure  driven  flows  of  propellants  can  give  rise  to 
flow  instabilities  associated  with  the  migration  of  the  polymeric  binder  in  the 
direction  of  the  pressure  gradient,  which  drives  the  flow  [21-24].  Such  demixing  of 
the  ingredients  of  the  formulation  of  highly  filled  suspensions  affects  the  quality  of 
the  degree  of  mixedness  of  the  suspension  and  hence  its  ultimate  properties.  It 
renders  the  continuous  processing  operation  unstable  and  also  poses  a  safety 
threat  in  processing  of  energetic  suspensions.  The  use  of  gear  pumps  or  sensors 
with  tight  converging  sections  also  become  prohibiting  for  propellants  and  other 
PEP’s. 

The  basic  driving  force  of  flow  instabilities  in  the  flow  of  propellants  is  the  axial 
migration  of  the  binder,  "filtering,"  superimposed  on  the  pressure-induced  bulk 
flow  of  the  energetic  suspension.  Upon  axial  migration  concentration  gradients 
develop  in  the  energetic  suspension  found  in  the  reservoir  under  pressure  and 
the  extrudate  emerging  from  the  die.  The  axial  migration  of  the  binder 
phenomenon  is  demonstrated  in  Figure  12  [24].  Here,  upon  the  pressure  induced 
flow  of  a  suspension  in  a  cylindrical  reservoir  connected  to  a  capillary  die,  the 
propellant  simulant  suspension  (used  in  ASRM  program  for  solid  rocket  motors) 
loses  some  of  its  binder  content  and  develops  a  concentration  gradient  in  the 
axial  flow  direction,  z,  between  the  plunger  which  drives  the  flow  and  the  die. 
Under  severe  unstable  flow  conditions  the  pressure  necessary  to  extrude  the 
propellant  oscillates  while  growing  unbounded  with  time.  The  time  periodic 
oscillations  in  extrusion  pressure  are  related  to  the  time  periodic  mechanism  of 
formation  of  a  mat  of  solids  at  the  capillary  or  die  and  its  break-up  which  occurs 
concomitant  with  the  filtration  of  the  binder  in  the  axial  direction  [22]. 

The  onset  of  flow  instabilities  associated  with  mat  formation  and  filtering  of  the 
binder  can  be  predicted  with  reasonable  accuracy  by  comparing  the  filtration 
velocity,  Vm  (average  velocity  of  the  binder  through  a  packed  bed  of  solids)  of  the 
binder  with  the  wall  slip  velocity,  Us,  under  the  generally  prevailing  plug  flow 
conditions  [22,24].  The  flow  generally  becomes  unstable  when  the  filtration 


Figure  12.  Concentration  gradient  development  in  a  density  matched 
suspension  (64.6%  solid  polymethyl  methacrylate  particles  in  polyether  oil 
based  binder)  upon  pressure  flow  in  a  cylindrical  reservoir  connected  to  a 
capillary  die  [24]. 


velocity  Vm  of  the  binder  becomes  greater  than  the  slip  velocity  of  the  propellant 
plug.  This  phenomenon  is  shown  in  Figure  13  for  the  ASRM  propellant  simulant 
containing  76.5  percent  vol.  solids.  Line  C-C'  represents  the  points  where  the  slip 
velocity,  Us,  is  equal  to  the  filtration  velocity  of  the  binder.  Theoretically,  extrusion 
at  apparent  shear  rates  found  to  the  left  of  line  C-C'  is  subject  to  the  time- 
periodic  filtering  of  the  binder  and  constitute  the  unstable  region.  The  range 
compares  favorably  with  the  experimentally  observed  thresholds  for  the  on-set  of 
flow  instabilities  represented  by  line  B-B’. 

Obviously  in  the  processing  of  propellants  the  filtration  of  the  matrix  and  thus  the 
unstable  region  should  be  avoided  [21-24].  Otherwise,  the  additional  loss  of  the 
binder  due  to  axial  migration  of  the  polymeric  binder  at  solid  loading  levels,  which 
are  very  close  to  the  maximum  packing  fraction  of  the  fuel  and  oxidizer  particles 
will  render  the  propellant  unprocessable.  The  unstable  region  in  converging  flows 
can  be  avoided  by  the  proper  selection  of  the  production  rate  and  the  geometry 
of  the  channel  through  which  the  suspension  is  flowing.  We  have  studied  some 
of  the  pertinent  factors  and  reported  on  their  effects  to  help  the  propellant 
industry  to  properly  select  their  geometries  and  operating  conditions  (14,  21-24). 


VI.  Air  Entrainment  Effects  in  Processing  of  Propellant  Suspensions 

Since  the  formulation  of  a  highly  filled  propellant  suspension  involves  a  volume 
fraction  of  solids  i.e.,  fuel  and  oxidizer  particles  which  approaches  the  maximum 
packing  fraction  of  the  solid  phase,  the  incorporation  of  even  small 
concentrations  of  air,  makes  a  significant  difference  in  the  rheological  behavior 
and  hence  the  processability  of  propellant  filled  suspensions  [25-29].  The  shear 
viscosity  of  the  energetic  suspension  decreases  and  wall  slip  velocity  values 
increase  with  the  incorporation  of  air  into  the  propellant  [28]. 

The  amount  of  air  entrained  into  the  energetic  suspension  increases  at  partially 
full  regions  in  comparison  to  suspension  samples  found  at  the  completely  full 
sections  of  processors  during  the  continuous  processing  of  propellants  and  their 
simulants  [26].  The  density  of  samples  collected  from  the  reversely  configured 
screw  sections  of  a  twin  screw  extruder  (which  necessitate  completely  full  mixing 
volume)  are  greater  than  those  of  the  samples  collected  from  the  partially  full 
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Figure  13.  Experimentally  (B-B1)  and  theoretically  (C-C) 
determined  unstable  flow  zone  in  capillary  flow  of  a 
suspension  with  76.5%  solids  [24]. 


sections  of  the  same  extruder  i.e.,  at  forwardly  configured  screw  sections  [26]. 
The  air  entrained  into  the  suspension  also  affects  the  development  of  the  slip 
layer  (binder  rich  region  found  adjacent  to  the  wall)  of  the  suspension  during  die 
flows  [30].  Under  moderate  die  pressures,  pockets  of  air  cling  to  the  wall 
momentarily  and  are  spread  out  to  be  later  dragged-on  by  the  bulk  of  the 
suspension  [30]. 

VII.  Mathematical  Modeling  of  Twin  Screw  Extrusion  of  Propellants 

Currently  twin  screw  extrusion  processing  is  used  in  various  processing 
applications  where  the  ingredients  of  a  propellant  formulation  are  interpersed  in 
each  other  (distributive  mixing),  the  physical  properties  of  one  or  more  of  the 
components  are  altered  (dispersive  mixing  such  as  a  decrease  in  particle  size), 
the  air  or  solvent  content  of  the  suspension  is  reduced  under  vacuum 
(devolatilization),  and  the  energetic  suspension  is  pressurized  (for  shaping  in  a 
die  or  mold);  all  within  the  confines  of  a  single  processor  with  properly  selected 
screw  and  barrel  elements  [31].  Such  versatile  continuous  processors  provide  a 
better  surface  to  volume  ratio  and  facilitate  real-time,  on-line  process  control  and 
easier  pollution  prevention  in  comparison  to  batch  mixers.  However,  due  to  the 
peculiar  nature  of  the  flow  and  deformation  behavior  of  propellants  with  very  high 
degree  of  fill  their  twin  screw  extrusion  processing  modeling  is  a  special 
challenge  [32-41].  For  example,  the  wall  slip  behavior  needs  to  be  incorporated 
into  the  mathematical  modeling  of  their  processing  behavior  [32,  33,  35,  38,  39]. 
Figure  14  shows  the  typical  effect  of  the  ratio  of  the  wall  slip  coefficient  at  the 
screw  surface  over  the  slip  coefficient  at  the  barrel  surface,  cp,  on  the  operating 
curve  of  a  single  screw  extruder  [32],  At  constant  mass  flow  rate  the 
pressurization  capability  decreases  with  increasing  slip  coefficient  ratio,  rp. 

Simulation  Techniques  Developed 

Flow  Equations 

The  governing  equations  of  the  flow  in  the  extrusion  flow  are  given  by  the 
general  equations  of  conservation  of  mass  and  momentum,  which  are  fully  elliptic 
and  three-dimensional.  For  regular  flighted  elements  simplifying  assumptions 
can  be  made  which  render  the  equations  parabolic  and,  therefore,  amenable  to 
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Figure  14.  Operating  curve  (pressurization  rate  versus  production  rate) 
of  a  single  screw  extruder  as  affected  by  the  ratio  of  wall  slip 
coefficients  at  screw  and  barrel  surfaces,  <p  [32]. 


stepwise  integration  in  the  z-direction  (primary  flow  direction)  from  prescribed 

upstream  conditions.  These  assumptions  are  as  follows: 

(i)  the  inertia  effect  is  negligible  because  of  the  highly  viscous  nature  of  the 
generally  polymeric  materials  used  in  extrusion  and  the  creeping  nature  of 
the  resulting  flow, 

(ii)  the  flow  is  steady, 

(iii)  the  screw  channel  can  be  unwound  in  a  helical  direction  of  the  screw, 

(iv)  the  velocity  derivatives  along  the  channel  are  small  compared  to  those  in 
the  transverse  directions  except  for  locations  that  are  close  to  the  entrance, 
and 

(v)  the  barrel  moves  relative  to  the  screw  at  an  angle  (i.e.  the  helix  angle)  to  the 
downchannel  direction  with  a  constant  velocity  given  by  the  screw  rotational 
speed  and  barrel  diameter. 

However,  for  kneading  disc  type  elements  the  flow  field  can  only  be  solved 

in  3-D. 

The  equations  of  conservation  of  mass  and  momentum  are: 


0 


(1) 


d  P 
dx 


a_ 

dx 


f 


2q 


V 


ax 


dy 


(2) 


d  P  d 
dy  ~  dy 


a  f 

fdu  du\ 

N) 

+ 

x  y 

dy  +  dx 

v  7  J 

dp 


m 


dz 


d_ 

dx 


f  3uzA 
^  3x 


dy 


5u,^i 


n 


dy 


(3a) 

(3b) 


where  the  asterisk  (*)  denoting  the  dimensionless  forms  of  the  variables  has 
been  suppressed  in  the  equations  for  convenience.  For  regular  flighted  screw 
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sections  the  parabolic  nature  of  these  differential  equations  is  preserved  through 
pmj  the  mean  viscous  pressure  which  has  been  defined  for  the  primary  flow 
direction,  z,  and  is  constant  for  any  x-y  plane.  The  pressure  gradient  in  the 
downchannel  direction  pressure  gradient,  dpm/dz,  which  is  a  function  of  the 
channel  direction  coordinate,  z,  for  temperature  dependent  shear  viscosity,  is 
determined  by  the  requirement  that  the  conservation  of  mass  constraint  on  the 
total  volumetric  flow  rate  is  satisfied,  i.e.: 

Q-j  uzdD  =  0  (4) 

D 

where  Q  is  set  by  the  specified  inlet  value  of  dpm/dz.  The  pressure  p  ,  is  allowed 
to  vary  in  the  transverse  directions  in  such  a  way  that  the  continuity  equation  (Eq. 
1)  is  satisfied.  For  regular  flighted  elements,  even  though  uz  varies  along  the 
downchannel  direction,  z,  it  may  be  assumed  in  Eq.  (1)  that  step  sizes  in  z 
direction  are  so  small  that  5uz/9z  is  negligible  compared  to  the  other  two  terms. 
Equations  (2a)  and  (2b)  can  now  be  handled  separately  from  the  channel 
direction  equation  (2c),  and  the  penalty  method  is  used  in  approximating  the 
pressure  p  . 


Boundary  Conditions  and  Implementation  of  Wall  Slip 

Applying  Bubnov-Galerkin's  method  to  Eqs.(2a)  and  (2b)  leads  to  the 
following  dimensionless  residual  equations: 


[  Vrc-TdD 
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n(n  *T)drs  = 
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(5) 


where  k  is  the  weighting  function  andj",  the  total  stress  tensor  is  as  modified  by 
the  parabolic  flow  assumption.  On  the  boundaries,  the  traction  n  -J  can  be 
decomposed  into  the  tangential  and  normal  components  to  obtain: 
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where  t  and  n  are  respectively  the  unit  tangent  and  the  unit  outward  normal 
vectors  to  the  solid  boundary.  Eqs.  (2a)  and  (2b)  are  required  to  be  satisfied  on 
the  boundary.  The  conditions  of  wall  slip  and  no  material  exchange  between  the 
fluid  and  the  boundary  in  dimensionless  form  can  be  expressed  as: 

t  -  (u  -  uS0|jd)  =  P*(nt:J  (7a) 


n  ■  (u  -  u  ..  .)  =  0 
solid' 


(7b) 


where  uc  is  the  solid  boundary  velocity,  u  is  the  fluid  velocity, 

p*  =  (m0con_1  /  Rs)p  is  the  dimensionless  Navier's  slip  coefficient,  i.e.,  p  =  0  gives 
no  slip  and  p  =  °o  is  perfect  slip.  Eq.  7a  is  the  Navier's  slip  condition  as  applicable 
to  the  two-dimensional  transverse  plane  [23]. 


Referring  now  to  Eq.  6  for  the  implementation  of  the  boundary  conditions, 
one  observes  that  while  the  slip  condition  easily  blends  with  it  by  providing  the 
tangential  component  of  the  traction,  the  normal  component  presents  a  difficulty 
as  it  is  not  recoverable  from  the  second  condition,  i.e.,  Eq.  7b,  which  gives  only 
the  normal  velocity.  In  essence,  the  numerical  difficulty  posed  by  a  slip  boundary 
condition  is  the  requirement  that  both  the  essential  (kinematic)  condition  (Eq.  7b) 
and  the  natural  (Navier's  slip)  condition  be  satisfied  at  every  boundary  nodes.  We 
have  transformed  the  essential  condition  of  normal  velocity  into  a  natural 
condition  [25]: 


nt  ;T  =  t-(u-UsgJid)  (8a) 

P* 

nn  :T  =  n  (U-Us0'id)  (8b) 

where  a  new  penalty  parameter,  Xs,  has  been  introduced.  As  Xs  -»  0,  Eq.  8b 

becomes  equivalent  to  Eq.  7b  since  the  normal  component  of  the  traction  must 
be  finite  at  the  boundary.  In  practice,  a  Xs  value  in  the  range  10‘10  -  10'18 

produces  a  normal  velocity  difference  of  the  same  order  of  magnitude.  Equations 
8a-b  are  now  substituted  into  Eq.  6  to  give: 


Jvjt  T  dD  -  |  .*•<"  “?olid*  tdTg-  j  )tn'*U^Usol'd*ndrs  =  0  (9) 

D  rs  P  rs 

This  method  is  generally  applicable  to  any  rheological  behavior,  geometrical 
shape  of  the  boundary,  and  orientation.  For  the  channel  direction  momentum 
equation,  the  normal  component  of  traction  does  not  exist-which  makes  the 

incorporation  of  wall  slip  by  the  Galerkin  method  straightforward.  The 
components  of  uso|jd  on  the  screw  root  surface  are  all  zero  while  on  the  barrel, 

the  solid  boundary  velocity  uso|id,  is  given  in  dimensionless  form  as: 

usolid  =  isino  +  kcoso  (10) 

where  0  is  the  helix  angle  of  the  screw  and,  i  and  k  are  unit  vectors  in  x-  and  z- 
directions  respectively. 

Constitutive  Model 

The  propellant  formulations  are  assumed  to  be  purely  viscous  fluids  which 
can  be  characterized  Theologically  by  the  generalized  Newtonian  fluid  constitutive 
relation: 


I  =  -V(7)y 


(11) 


where  r)  is  the  shear  viscosity  material  function,  y  is  the  rate-of-deformation 
tensor  and  y  is  the  deformation  rate.  The  dependence  of  the  shear  viscosity 
material  function  on  the  deformation  rate  and  temperature  follows  the  modified 
Herschel-Bulkley  model: 
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where  n  is  a  material  parameter  which  governs  the  sensitivity  of  the  fluid  to 
deformation  rate,  xy  is  the  apparent  yield  stress,  nb  is  the  stress  growth 
exponent,  c'is  the  temperature  coefficient  of  viscosity,  with  To,  the  entrance 


temperature  as  the  reference  temperature.  The  magnitude  of  the  rate-of- 
deformation  tensor  |y|  is  defined  for  our  case,  by: 


“  ^ux,x  +  ^uy,y  +  (ux,y  +  uy,*f  +  uz,x  +  uz,y 


(13) 


where  comma  indicates  differentiation.  If  we  introduce  the  following 
dimensionless  variables: 


|y|  =  <o|y|*,Ty  =  mo0)nTy,nb  =  nb  /  co 


(14) 


then  Eq.  12,  in  dimensionless  form,  becomes: 
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Various  simplifications  of  the  modified  Herschel-Bulkley  model  include  the 

*  * 

Ostwald-de  Waele  power-law  fluid  (xy  =0),  the  Newtonian  fluid  (n=1,  ty  =0), 

the  modified  Bingham  plastic  (Papanastasiou  [26],  (n=1)  and  the  classical 

* 

Bingham  plastic  fluid  above  the  yield  point  (n=1,  nb=°°).  For  any  given 
propellant  formulation  the  parameters  of  the  constitutive  equation  need  to  be 
determined  using  the  viscometric  flows  discussed  earlier.  The  Navier's  wall  slip 
coefficient  p*  also  needs  to  be  characterized  for  a  given  propellant  from 
viscometric  experimental  techniques  which  we  have  developed  (8,  9,  12,  13,  16). 


Energy  Equation  and  Boundary  Conditions 

By  an  order-of-magnitude  analysis  of  the  energy  equation,  it  can  be 

demonstrated  that  in  the  extrusion  flow  the  importance  of  the  transverse 
convection  terms  (pCpux3T/ax,  pCpuy3T/ay)  relative  to  heat  conduction  in  y- 

direction  scales  as  Peu*  H2/W  where  Pe  =pC  coR?/k  is  the  Peclet  number  and  H 

X  P  o 

and  W  are  respectively,  the  dimensionless  maximum  depth  and  width  of  the 

screw  channel.  Polymeric  binders  and  propellants  generally  have  low  thermal 
conductivity  and,  moderately  high  values  of  p  and  Cp,  hence  exhibit  high  Peclet 

number.  Therefore,  the  transverse  convection  term  will  be  negligible  only  for 


small  H  and  large  W  (i.e.  low  aspect  ratio  geometries)  and  also  for  small  helix 
angle.  If  the  transverse  convection  term  is  compared  with  the  downchannel 
convection  term,  the  transverse  convection  term  will  be  negligible  if  the  ratio 
uxz* /uz\N  is  small,  where  u*  and  are  dimensionless.  For  the  regular  flighted 

element,  from  the  continuity  equation  and  the  typical  choice  of  the  step  size, 
u*  /z*  is  generally  small.  Hence,  the  transverse  convection  terms  can  be  ignored 

only  if  W  is  very  large  or  the  helix  angle  is  very  small — the  same  conditions  as 
obtained  earlier.  Since  these  conditions  are  rarely  encountered  in  the  design  and 
operation  of  realistic  extruders,  the  transverse  convection  terms  should  be 
included  in  the  energy  equation  for  the  analysis  of  extrusion  flows. 


By  a  similar  analysis,  the  quantity  Peu*  z  governs  the  relative  importance 

of  the  downchannel  convection  to  heat  conduction  in  the  downchannel  direction. 
Therefore,  except  for  low  volume  flow  rates  and  the  region  very  close  to  the 
entrance,  one  can  neglect  heat  conduction  in  the  downchannel  direction  since 
the  Peclet  number  is  generally  high.  With  these  assumptions  and  the  assumption 
of  steady  state  conditions,  using  the  dimensionless  temperature  0  =  (T  -  T0)  /  (Tb 
-  T0),  the  energy  equation  in  dimensionless  form  becomes: 
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where  G=m0con  1  (coRs )2/k(Tb-  TQ)  is  the  Griffith  number  (generally  referred  to 

as  the  Brinkman  number  for  the  Newtonian  fluid).  The  asterisk  (*)  denoting  the 
dimensionless  forms  has  been  suppressed  for  the  velocity  components  and  co¬ 
ordinates.  The  Griffith  number  indicates  the  relative  importance  of  viscous 
dissipation  effects  to  heat  conduction. 

The  fluid  enters  the  screw  channel  at  a  temperature  TQ.  The  barrel 
temperature  is  designated  as  Tb  and  is  different  from  the  entrance  temperature. 

Variations  in  Tb  along  the  channel  direction  are  allowed.  Although  our  model  can 
be  used  in  conjunction  with  any  type  of  thermal  boundary  condition,  we  have 
found  that  constant  barrel  temperature  and  adiabatic  screw  surface  conditions 
are  proper  boundary  conditions.  Thus,  in  dimensionless  form,  the  initial  and 
boundary  conditions  become: 
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at  z  =  0,  0  =  0 
on  barrel  surface  0  =  1 

on  screw  surface  30  /  3n  =  0  (17a-c) 

where  n  is  the  unit  vector  normal  to  screw  surface.  Equation  (16)  was  initially 
solved  with  the  above  specified  boundary  conditions  using  the  Bubnov-Galerkin 
method  and  the  removal  of  oscillations  in  the  temperature  solution  was  initially 
attempted  by  a  simple  smoothening  procedure.  However,  later  we  modified  our 
algorithms  to  implement  the  recently  developed  SUPG  method  (which  has  rarely 
been  implemented  in  polymer  extrusion  simulation)  by  introducing  weighting 
functions  of  the  form: 

Wf  =  Wf  +  wf  (18) 

where  Wf  is  a  continuous  function  of  space  and  wf  is  a  discontinuous  perturbation 
of  Wf.  The  weighted  residual  form  of  the  energy  equation  is  then  given  by: 

dD 

-Jr  Wf(V0-n)drs|dz  =  O  (19) 

The  shape  functions  are  chosen  to  satisfy  the  boundary  conditions  on  the  barrel 
and  must  therefore  vanish  there. 

Several  approaches  have  been  suggested  in  the  literature  for  the  definition 
of  the  perturbation  function  for  SUPG.  We  have  followed  the  method  outlined  by 
Hughes  et  al.  In  this  analysis,  an  element  in  the  arbitrary  domain  x,  y,  and  z  can 
be  mapped  into  an  isoparametric  element  of  unit  dimension  in  each  of  the 
transformed  serendipity  coordinates  %,  and  An  element  mesh  parameter  was 
then  defined  as 
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where  Hull  is  the  Euclidean  norm  of  the  velocities  and  b  is  given  by 
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Following  Hughes  et  al.  the  perturbation  function  was  taken  as: 
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where, 
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NUMERICAL  SOLUTION 

Flow  Equations 

For  the  finite  element  solution  of  the  momentum  equations,  the  domain  is 
discretized  by  bilinear  elements  in  x  and  y  and  trial  functions  are  chosen  such 
that: 

Nod 

xe  =  X  xfnf  (x,y)  xf  ,nf  e  De  (24  a,  b) 

i  =  1 

where  De  is  an  element  domain,  Nod  is  the  number  of  nodes  per  element  and  n ® 

is  the  trial  function  with  xe  representing  the  velocity  vector.  The  discretized 
equations  take  the  matrix  form: 


A  x  =  f  (25) 

where  A  is  the  stiffness  matrix  and  f  is  a  force  vector  resulting  from  specified 
nodal  values  or  source  term  and  the  nodal  unknowns  are  recovered  from  the 
vector  x. 


Energy  Equation 

The  shape  functions  for  the  temperature  are  similar  to  those  used  for 
velocity  in  that  they  are  linear  in  the  space  coordinates.  The  weighting  function 
Wf  has  a  trilinear  dependence  on  the  x-,  y-,  and  z-  coordinates  for  regular  flighted 
screw  sections.  Since  the  thermal  problem  is  parabolic,  the  integration  in  z  in  Eq. 
16  can  be  carried  out  independent  of  the  other  two  directions  producing  a  part 
with  known  values  of  9  (source  term)  and  the  other  containing  values  of  9  being 
sought.  At  this  point,  the  solution  procedure  adopted  for  the  velocity  equations 
can  then  be  applied. 

The  element  mesh  parameter  If  and  the  cell  Peclet  number  ye  in  Eqs.  20 
and  23  were  evaluated  based  on  the  velocity  at  the  centroid  of  the  element.  Our 
numerical  experimentation  revealed  that  the  value  of  a  is  dependent  on  such 
factors  as  the  element  skewness  and  the  flow  and  heat  transfer  parameters.We 
have  recommended  that  a  be  adjusted  until  the  solution  is  oscillation  free.  For  the 
range  of  parameters  considered  in  various  case  studies,  the  values  of  a  were 
between  0.5  and  3.0  and  they  furnished  satisfactory  results. 

Solution  Procedure 

For  regular  flighted  screw  sections,  at  the  inlet  of  the  extrusion  channel, 
the  temperature  distribution  is  known,  hence  only  the  x-,  y-,  and  z-  velocity 
components  are  computed,  based  on  a  specified  value  of  the  inlet  pressure 
gradient  dpm*/dz*.  The  Picard  iteration  is  then  utilized  for  the  shear  rate 
dependent  shear  viscosity.  After  convergence,  the  volumetric  flow  rate  is 
calculated.  We  proceed  down  the  channel  and  as  a  first  step,  solve  for  the 
temperature  distribution.  With  a  viscosity  distribution  based  on  this  temperature 
profile  and  the  value  of  dpm7dz*  at  the  previous  channel  location,  for  the  second 
step,  we  solve  for  the  velocity  components  and  calculate  the  volumetric  flow  rate. 
A  new  value  of  dpm7dz*  is  obtained  using  the  secant  method  with  the  residual 
taken  as  the  difference  between  the  calculated  volumetric  flow  rate  and  that 
obtained  at  the  inlet.  This  completes  the  second  step  which  is  repeated  until 
convergence  on  account  of  shear  rate  dependent  shear  viscosity. 

Subsequently,  we  solve  for  a  new  temperature  distribution  and  go  back  to 
the  second  step  since  the  shear  viscosity  material  function  is  temperature 
dependent.  This  iterative  procedure  is  continued  until  the  volumetric  flow  rate 
constraint  (Eq.  4)  is  satisfied  and,  simultaneously,  the  Euclidean  norms  of  the 


difference  between  successive  iterates  of  u  and  0  are  small,  usually  <  10'3. 
Another  step  size  is  taken  with  the  values  of  the  variables  at  the  previous  channel 
location  serving  as  the  initial  guesses.  Convergence  on  the  inner  iteration  for 
shear  rate  dependent  shear  viscosity  is  usually  achieved  in  4  to  8  iterations.  On 
the  other  hand,  the  outer  Picard  iteration  on  account  of  the  temperature 
coefficient  of  viscosity  may  require  up  to  10  iterations.  The  number  of  iterations 
depends  on  the  channel  location  and  the  value  of  the  temperature  coefficient  of 
viscosity.  When  the  backflow  is  significant  but  not  severe,  the  SUPG  method 
provides  reliable  and  accurate  results  when  compared  with  the  upwind  marching 
scheme  of  finite  difference,  but  the  required  number  of  iterations  goes  up  slightly. 
The  allowable  step  size,  Az,  is  dictated  by  the  stability  of  the  discretized 
temperature  equation  which  is  dependent  on  the  thermal  and  flow  parameters. 
The  dimensionless  step  size  used  in  this  study  is  1.0. 

For  the  kneading  disc  sections  of  the  twin  screw  extruder  the  full  three- 
dimensional  conservation  equations  need  to  be  solved  and  the  procedures  used 
are  more  complicated.  They  are  provided  in  detail  in  some  of  our  publications 
(34,  37,  40,  41). 

The  typical  role  played  by  the  viscoplasticity  of  highly  filled  materials  on  their 
continuous  processability  is  demonstrated  in  Figures  15,  16.  Here  the 
suspension  is  pressurized  in  a  co-rotating  twin  screw  extruder  [34,  37,  40-42]. 
The  velocity  distributions  given  in  Figure  15  show  that  the  flow  is  segregated  into 
a  high  deformation  rate  region  adjacent  to  the  barrel  surface  and  a  low 
deformation  rate  zone  away  from  the  barrel.  The  consequence  of  this  flow  field  is 
the  formation  of  a  hot  spot  adjacent  to  the  barrel  surface  as  shown  in  Figure  16, 
which  needs  to  be  avoided  by  proper  selection  of  geometry  and  processing 
conditions. 

An  additional  difficult  aspect  of  continuous  processing  arises  from  the  dispersive 
mixing  associated  with  the  exertion  of  high  stress  levels  on  small  fractions  of  the 
suspension  by  forcing  repeatedly  through  the  small  gaps  found  between  screws 
or  the  screw  and  the  barrel  elements  during  continuous  processing  [34,  40].  This 
is  illustrated  in  Figure  17  where  high  stress  magnitudes  are  seen  to  be  generated 
at  such  gaps  during  the  processing  of  a  viscoplastic  suspension  in  a  co-rotating 
twin  screw  extruder.  The  physical  properties,  including  particle  size  distribution, 
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Figure  17.  Distribution  of  the  shear  stress  in  a  viscoplastic  fluid  being 
processed  in  the  kneading  disc  region  of  a  co-rotating  twin  screw 
extruder. 


of  the  suspension  experiencing  the  high  stress  region  can  become  significantly 
different  than  the  rest  of  the  suspension. 

Degree  of  Mixedness 

Some  of  our  experimental  and  numerical  simulation  techniques  which  were 
developed  earlier  with  DARPA  funding  were  subsequently  applied  to  the 
investigation  of  propellant  processing  in  co-rotating  twin  screw  extrusion  process 
using  earlier  funding  provided  by  DARPA.  Besides  simulation  technologies  [36, 
37,  40,  41],  we  have  developed  two  experimental  techniques  which  are  both  x- 
ray  analysis  based  to  characterize  the  volume  ratio  of  various  ingredients  as  a 
function  of  the  scale  of  examination  and  location  [43].  The  techniques  are 
illustrated  in  Figure  18  where  the  distributive  mixing  of  a  minor  component 
indicated  by  a  white  blob  in  the  geometry  of  two  co-rotating  kneading  discs 
confined  in  a  barrel  is  shown.  On  the  other  hand,  we  have  also  developed 
simulation  techniques  for  predicting  the  degree  of  mixedness  development  in 
propellants  during  extrusion  flow. 

Figure  19  illustrates  an  x-ray  based  experimental  technique,  which  we  have 
developed  [43].  The  variation  of  the  volume  ratios  of  the  binder,  cellulose  acetate 
butyrate,  (CAB)  over  the  solid  content,  i.e.,  RDX  and  HMX  powders  at  various 
locations,  are  shown.  Samples  were  either  processed  in  a  batch  mixer  and  then 
ram  extruded  or  processed  and  extruded  in  a  twin  screw  extruder.  Figure  20 
shows  that  the  volume  fraction  of  the  binder,  CAB,  is  seen  to  vary  considerably 
around  the  mean  for  the  batch  mixed  sample  versus  the  narrow  distribution 
around  the  mean  obtained  with  the  continuously  processed  gun  propellant 
sample.  The  difference  becomes  more  significant  at  a  smaller  scale  of 
examination  i.e.,  at  1  mm^-  Predictions  assume  that  the  minor  component  is 
non-diffusive  with  no  interfacial  tension  existing  between  the  minor  component 
and  the  remainder  of  the  suspension  [36].  As  also  experimentally  observed,  the 
distribution  of  the  minor  component  in  the  remainder  of  the  suspension  becomes 
uniform  with  time  [44],  We  utilize  the  variance  of  the  ratio  of  the  minor  to  the 
major  phases  as  a  quantitative  measure  of  degree  of  mixedness.  The  challenge 
is  to  determine  proper  locations  for  minor  components  to  be  input  into  the 
processor  and  selection  of  appropriate  geometry  and  operating  condition  for 
efficient  mixing  without  detoriating  various  properties  of  the  suspension. 


Figure  1 8.  Distributive  mixing  of  a  non-diffusive  passive  tracer  blob  as  a  function 
of  time,  t,  i.e.,  in  periods  of  rotation,  to,  2-D  flow  between  two  kneading  disks  [36]: 
a)  t  =  0;  b)  t  =  2  periods;  c)  t  =  5  periods;  d)  t  =  10  periods;  e)  t  =  20  periods. 
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Figure  19.  Extruded  profile  and  sampling 
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Figure  20.  Quantitative  analysis  of  degree  of  mixedness  using  the  variance 
of  the  concentration  of  the  binder,  CAB,  at  two  scales  of  examination  for 
batch  and  continuously  processed  suspensions. 


The  mixing  index  values  range  from  zero  for  a  completely  segregated  sample 
(unmixed)  to  one  for  a  perfect  mixture  with  zero  variance  around  the  mean.  The 
mixing  index  values  for  the  batch  mixed  sample  of  Figure  20  is  0.85  and  the  twin 
screw  processed  sample  is  0.92  at  a  scale  of  examination  of  1  mm2,  indicating 
the  superior  distributive  mixing  capability  of  the  twin  screw  extrusion  process. 
This  analysis  technique  and  the  simulation  techniques  for  the  first  time  allow  the 
quantitative  characterization  of  the  degree  of  mixedness  of  suspension  samples 
and  facilitate  the  determination  of  proper  mixing  geometries  and  conditions  for 
propellants  in  twin  screw  extrusion  process. 

Figure  21  shows  the  distribution  of  the  tracer  particles  (numerical  tracer  i.e., 
absence  of  interfacial  effects)  initially  forming  a  coherent  blob  located  at  the  entry 
point  of  a  kneading  disc  section  of  the  co-rotating  twin  screw  extruder.  The  rapid 
redistribution  of  the  tracer  particles  with  increasing  residence  time  in  the 
kneading  section  was  characterized  by  the  first  solution  of  the  time-dependent 
velocity  distributions  in  the  kneading  disc  section  followed  by  determination  of 
particle  paths. 

Figure  22  shows  the  typical  residence  time  distribution  plots  for  various  kneading 
disc  configuration  in  the  co-rotating  twin  screw  extrusion  process.  The  effect  of 
decreasing  of  the  stagger  angle  is  to  increase  the  breadth  of  the  residence  time 
distribution.  The  short  residence  time  values  at  30°  stagger  are  indicative  of  the 
"pipeline"  flows  (which  were  also  experimentally  observed  for  30°)  and  the  long 
residence  times  support  the  formation  of  relatively  stagnant  regions  in  the 
extruder. 

VIII.  Free  Surface  Flows 

A  number  of  sections  of  the  twin  screw  extruder  operate  on  the  basis  of  partially- 
full  mode.  This  is  typically  shown  in  Fig.  23  where  the  degree  of  fill  changes  with 
the  geometry  employed  (increase  with  reversely  configured  sections  and 
decrease  as  the  stagger  angle  decreases  for  forwardly  configured  screw 
sections).  For  the  regular  flighted  screw  sections,  the  partially  full  section  only 
serves  as  a  conveying  geometry  with  little/no  distributive  or  dispersive  mixing  of 
the  ingredients.  The  degree  of  fill  in  such  sections  and  the  residence  time  can  be 


Figure  21.  Mixing  of  a  blob  for  kneading  discs  staggered  at 
90°  netural:  (a)  t  =  0,  (b)  t  =  T0,  (c)  =  5T0,  and  (d)  t  =  20To 

neutral. 


Figure  22.  The  comparison  of  the  residence  time 
distribution  functions  generated  by  kneading  disc 
screw  sections  staggered  at  30°,  45°,  60°  forward, 
and  90°  neutral. 
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Figure  23.  Degree  of  fill  distribution  with  various  kneading  disc  stagger  angles  at 
70°C  and  40  rpm. 


determined  easily  since  the  volumetric  flow  rate  for  pure-drag  flow  is  not  a 
function  of  the  rheological  behavior  of  the  suspension  nor  the  wall  slip  condition. 

On  the  other  hand,  the  partially-full  sections  in  the  kneading  disc  sections  cannot 
be  modeled  with  the  present  state  of  knowledge.  In  such  sections,  the  motion  is 
governed  by  chaotic  dynamics  with  positive  Lyapunov  exponents  and 
characterized  by  absolute  dependence  on  the  initial  condition. 

The  regular  flighted  screw  sections  found  in  the  devolatilization  section  are  where 
the  mass  transfer  area  i.e.,  interfacial  area  between  air/vacuum  and  the 
suspension  is  important.  We  have  applied  the  "method  of  spines"  technique 
developed  by  Prof.  Scriven  and  co-workers  to  determine  the  shape  of  the  free 
surface  as  a  function  of  surface  tension,  contact  angles  and  wall  slip  coefficient  of 
the  suspension  on  the  extruder  surfaces.  Typical  results  for  x-and  z-velocity 
(transverse  and  down-channel  directions)  are  shown  in  Figure  24,  where  the  free 
surface  is  indicated.  The  determination  capability  of  the  free  surface  shape, 
however,  does  not  immediately  translate  to  a  devolatilization  model,  since  the 
nucleation  and  growth  of  bubbles  and  their  collapse  during  devolatilization  are 
phenomena  which  are  hitherto  not  mathematically  modeled.  We  have  attempted 
such  a  model  in  conjunction  with  development  of  pore  formation  in  conversion  of 
polymeric  precursors  to  ceramics  (Ref'ns  46),  but  had  no  opportunity  to  adapt 
and  carry-over  to  devolatilization  process. 

The  ramifications  of  all  these  findings  and  capabilities  to  the  BMDO  program  in 
terms  of  reducing  the  safety  risk  in  continuous  processing  and  optimization  of  the 
process  towards  "total  quality  management"  through  the  understanding 
generated  and  the  numerical/experimental  tools  developed  are  discussed  next. 


IX.  Conclusions  and  Ramifications 

1 .  The  flow  and  deformation  behavior  of  propellants  is  affected  by  the  formation 
of  a  binder-rich  region  adjacent  to  the  wall  i.e.,  wall  slip  region.  Generally, 
the  energetic  suspension  flows  like  a  plug  (solid-like),  lubricated  at  the  wall 
with  the  slip  region  adjacent  to  the  wall. 


Mathematical  modeling  of  2-D  free  surface  flows 
using  “method  of  spines.” 
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SIT’s  mathematical  model: 

Determines  free  surface  shape  in  2-D  flows 

Predicts  the  velocity,  stress  and  temperature  distributions  in  partially  full  channel 
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2.  The  solid-like  behavior  of  propellants  in  various  flow  regimes  renders  their 
processing  in  twin  screw  extrusion  processing  i.e.,  mixing  and  pressurization 
and,  shaping  in  dies  and  molds  a  challenge. 

3.  On-line  analysis  techniques  attempting  to  probe  various  chemical  or  physical 
aspects  of  propellants  should  be  able  to  penetrate  beyond  the  wall  slip  layer 
to  characterize  the  flowing  suspension. 

4.  The  thickness  of  the  deforming  slip  layer  changes  as  a  function  of  binder  and 
solid  phase  characteristics  of  the  propellant. 

5.  Plug-like  flow  of  the  propellant  occurs  at  low  wall  shear  stress  values  for 
shear  thinning  propellant  formulations  with  multimodal  size  particles  and  at 
high  wall  shear  stress  values  for  shear  thickening  propellant,  with  narrow 
particle  size  distributions. 

6.  Formation  of  solid-like  plugs  with  lubrication  layers  at  the  wall  renders  it 
difficult  to  achieve  adequate  mixing  and  processing  of  propellants  requiring 
specialized  mixing  and  dispersion  elements  in  twin  screw  extrusion. 

7.  The  microstructure  of  the  propellant  suspension  evolves  continuously  as  a 
function  of  flow  and  deformation  history.  Microstructure  development 
encompasses  migration  of  particles  based  on  shear  rate  and  concentration 
gradients,  interlocking  of  particles  to  form  solid-like  plugs,  encapsulation  of 
binder  locally  upon  interlocking,  formation  of  binder-rich  layers  at  high  shear 
rate  regions,  entrainment  of  air/solvent  gases,  and  lubrication  of  surfaces 
with  air/gases. 

8.  Especially  with  propellant  suspensions  containing  low  viscosity  binders,  the 
binder  migrates  in  the  direction  of  the  overall  pressure  gradient,  with  such 
filtration  superimposed  on  the  bulk  flow  of  the  propellant  suspension.  This 
migration  results  in  the  demixing  of  the  propellant  upon  pressurization  in  die 
and  mold  flows  and  leads  to  cyclic,  mat-formation  based  flow  instabilities. 

9.  The  factors  involved  in  binder  migration  and  development  of  flow  instabilities 
are  understood,  with  major  effects  arising  from  degree  of  channel 


convergence,  particle  diameter,  and  binder  viscosity.  Binder  migration  during 
propellant  processing  should  be  minimized  to  eliminate  the  formation  of 
stagnant  solid  regions  in  processors,  particle  attrition,  quality  problems  and 
increasing  safety  risks  during  continuous  processing. 

10.  Attempts  to  eliminate  wall  slip  during  deformation  by  roughening  the  walls  of 
processors  generate  fracture  surfaces  (discontinuities  of  velocity)  in 
propellants. 

11.  The  development  of  wall  slip  behavior  of  propellants  during  deformation  is 
time  and  temperature  dependent. 

12.  The  peculiarities  of  the  flow  and  deformation  behavior  of  propellant 
suspensions  with  their  high  degree  of  fill  require  specialized  numerical 
analysis  techniques  for  mathematical  modeling  of  their  processability  in  twin 
screw  extrusion  process  which  especially  incorporate  the  wall  slip  behavior. 

13.  The  wall  slip  behavior  of  propellant  suspensions  can  be  characterized  in 
viscometric  flows  and  used  for  first-order  analysis  of  2-D  and  3-D  processing 
flows  as  an  interfacial  dynamic  boundary  condition. 

14.  The  rheological  behavior  of  the  propellant  suspensions  is  especially  affected 
by  the  presence  of  a  gaseous  phase  i.e.,  air  or  solvents  and  is  altered 
significantly  upon  their  removal.  The  gas  entrainment  characteristics  of  the 
propellant  suspension  differs  significantly  at  different  locations  in  a 
continuous  processor  especially  varying  between  the  partially-full  and 
completely  full  sections  of  the  twin  screw  extruder. 

15.  The  wall  slip  behavior  and  hence  processability  of  propellant  suspensions 
can  be  affected  to  some  extent  by  the  choice  of  materials  of  construction  with 
different  material  use  on  moving  and  stationary  surfaces  of  the  twin  screw 
extruder,  being  the  most  effective. 

16.  The  microstructure  of  propellant  suspensions  which  controls  their  flow  and 
deformation  behavior  as  well  as  their  various  ultimate  properties  like  burning 
rate  and  mechanical  properties  are  affected  foremost  by  the  degree  of 


interspersing  of  various  ingredients  of  the  formulation  and  the  physical 
properties  of  the  ingredients  (such  as  particle  size  of  solids  which  can 
change  during  processing). 

17.  On-line  rheometers  for  highly  filled  materials  such  as  propellants  should  be 
able  to  probe  at  a  number  of  deformation  rates  without  significantly  altering 
the  microstructure  of  the  suspension. 

18.  The  axial  migration  of  binder  in  converging  flow  generally  prevents  the  use  of 
gear  pumps  and  sensors  where  suspension  is  dragged  into  small  clearances 
for  highly  filled  materials. 

19.  Small  changes  in  characteristics  of  raw  ingredients  such  as  particle  size  and 
shape  of  the  solid  phase  of  propellant  suspensions  where  the  intended  solid 
volume  concentration  approaches  the  maximum  packing  fraction  may  result 
in  suspensions  which  are  not  processable  (for  example  arising  from  the 
segregation  of  sold  ingredients  according  to  size  or  density  in  bins  and 
hoppers). 

20.  X-ray  based  techniques  developed  earlier  could  be  also  applied  to 
characterize  the  degree  of  mixedness  of  suspensions  processed  in  co¬ 
rotating  twin  screw  extrusion  process. 

21.  The  processability  of  concentrated  suspensions  could  be  simulated 
assuming  a  continuum  by  solving  the  conservation  equations  using  2-D  and 
3-D  Finite  Element  Method. 

22.  The  mixing  effectiveness  of  twin  screw  extrusion  geometries  and  operating 
conditions  could  be  quantitatively  assessed  using  tools  of  dynamics 
employing  non-diffusive  numerical  tracers. 

23.  The  generic  understanding  gained  by  applying  multi-disciplinary  studies  on 
propellants  and  their  simulants  is  also  relevant  to  other  materials 
encountered  in  multiple  industries  relevant  to  Navy. 


24.  New  technologies  can  be  developed  or  existing  technologies  improved  by 
understanding  the  generic  nature  of  propellant-like  materials  which  are  filled 
very  close  to  their  maximum  packing  fraction. 
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